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We apply a singularity resolution technique utilized in loop quantum gravity to the polymer 
representation of quantum mechanics on R with the singular — potential. On an equispaced 
lattice, the resulting eigenvalue problem is identical to a finite difference approximation of the 
Schrodinger equation. We find numerically that the antisymmetric sector has an energy spectrum 
that converges to the usual Coulomb spectrum as the lattice spacing is reduced. For the symmetric 
sector, in contrast, the effect of the lattice spacing is similar to that of a continuum self-adjointness 
boundary condition at a; = 0, and its effect on the ground state is significant even if the spacing is 
much below the Bohr radius. Boundary conditions at the singularity thus have a significant effect 
on the polymer quantization spectrum even after the singularity has been regularized. 



I. INTRODUCTION 



It is expected that a viable quantum theory of gravity will have to say something about what happens to the 
I , curvature singularities in classical solutions of general relativity. This would involve making concrete intuitive ideas 
^jrj' about the role to be played by the Planck length as a fundamental discreteness scale. A possible guide concerning how 
—I, to do this is the quantum resolution of the Coulomb potential in quantum mechanics, where the basic result is that the 
^ ' expectation value is finite in all eigenstates of the Hamiltonian. This kinematic result together with dynamical 



Coulomb scattering may be taken to constitute quantum singularity avoidance associated with the classically singular 
' —1/r potential. 

[ — , In the Hamiltonian approach to quantum gravity, the problem of quantization from a mathematical point of view 
04 ' is to find a suitable representation of an algebra of functions of position and momentum as operators on a Hilbert 
I space. The Wheeler-DeWitt approach utilizes a functional Schrodinger representation where the basic variables 
are the Arnowitt-Deser-Misner (ADM) variables [J, ■ The loop quantum gravity (LQG) approach uses a Poisson 
' algebra of functions of a connection and triad, based on loops and surfaces, to build a quantum theory [s^]. The basic 
; variables quantized are the holonomy of the connection along a curve, and the integral of the (densitized) triad over 
■ a surface. The latter approach naturally leads to a non-separable kinematical Hilbert space associated with graphs 
embedded in a spatial manifold Because of the association of kinematic states with graphs, there is an intrinsic 
spatial discreteness built into the quantum theory which is not a priori present in the Schrodinger approach. Field 
excitations of the basic operators are probed on graphs rather than at points. 



. An approach similar to that of loop quantum gravity can be employed for the quantum mechanics of a particle 
■ " moving in a potential in one or more dimensions. Here graphs are replaced by lattices of spatial points (not necessarily 
equispaced), and the basic observables realized in the quantum theory are the configuration and translation operators. 
This is of course natural since generators of infinitesimal translations cannot be represented on a spatial lattice. For 
a given equispaced lattice, the position eigenstates are normalizable, and the Hilbert space is obviously separable. 
However, the Hilbert space that incorporates all possible lattices (equispaced or not) is non-separable. The quantum 
theory that utilizes this space has been referred to as polymer quantization ^] . 

From a mathematical viewpoint, conventional Schrodinger quantization and polymer quantization are inequivalent. 
In the former, wave functions are square-integrable functions on M.^, with position and momentum operators acting 
as multiplication and differential operators, respectively. A remarkable fact is that this is quantization is unique up to 
unitary equivalence provided that the configuration space of the system is topologically M'^ and that the representation 
of the Weyl-Heisenberg algebra of exponentiated position and momentum operators is weakly continuous It is 



'Electronic address; lhusain@math.unb.cal 

t Electronic address: j orma. louko@nottingham. ac.ukj 

t Electronic address: owinkler@perimeterinstitute. ca| 



2 



known, however, that if either of these assumptions is abandoned, there are infinitely many inequivalent representa- 
tions. A simple example is a system where the configuration space is the torus rather than R"^. Polymer quantization 
provides another example. Its non-separable Hilbert space may be viewed as the inductive limit of the separable 
Hilbert spaces associated with quantum mechanics on all possible lattices, including the one-parameter family of 
equispaced ones. 

From a practical viewpoint, since only a finite number of calculations are possible, the full polymer Hilbert space 
can never be utilized. Rather, only a separable subspace is computationally useful. In this case, however, polymer 
quantization appears to reduce to the finitely-differenced Schrodinger equation on a lattice. We elaborate on this below 
by observing that various finite difference schemes for differential equations may be rewritten using configuration and 
translation operators on a lattice. 

In the representation used in LQG, there is a way to write inverse triad operators using certain classical Poisson 
bracket identities due to Thiemann These identities are in fact much more general than the context in which 
they first arose; similar ones may be written for any theory on a lattice. They may be used to write inverse scale 
factor operators that are bounded on kinematical states in quantum cosmology in both the connection-triad and 
ADM variables @, 0] j as well as curvature operators for a field theoretic formalism for gravitational collapse Q . The 
boundedness property of such operators has been used in discussions of singularity avoidance in quantum gravity at 
the kinematical level [1, [13] . 

In this paper we apply these ideas of singularity avoidance in quantum gravity to polymer representation quantum 
mechanics on K with the singular — l/|a;| potential. For x > 0, and with an appropriate boundary condition at 
X — 0, this may be thought of as the spherically symmetric sector of the Coulomb problem on R^. We address two 
main questions. First, we show that in any lattice context, the LQG singularity resolution technique is equivalent 
to replacing a singular derivative by a non-singular finite difference scheme. Second, we show that even after the 
singularity in the potential has been resolved, the spectra in the symmetric and antisymmetric sectors are significantly 
different. The latter is close to the usual Coulomb spectrum when the lattice spacing is much below the Bohr radius, 
but in the former the lattice spacing plays a role similar to a continuum self-adjointness boundary condition a.t x — 0, 
and the effect on the ground state is significant even when the spacing is much below the Bohr radius. We conclude 
that boundary conditions at the singularity have a significant effect on the polymer quantization spectrum even after 
the singularity itself has been regularized. 

The rest of the paper is as follows: In section |TT] we recall the basic structure of polymer quantization on M. In 
section imi we specialize to the potential — l/|a;|, introducing the lattice regularization of the potential and showing that 
the boundary conditions of the radial Coulomb problem can be implemented by the restriction to the antisymmetric 
sector. Our numerical results for the spectrum are given in section HVl The symmetric sector is analyzed in section fVl 
Section IVII summarizes the results and discusses their implications and limitations for the problem of singularity 
resolution in quantum gravity. 



In this section we briefly describe how a mechanical system is quantized in the polymer representation [j| . To keep 
the notation simple, and because the radial Coulomb problem we discuss later is a one-dimensional system, we will 
consider a particle on the real line. The generalization to n particles in MP is straightforward. 

Recall that the Hilbert space for the Schrodinger quantization of a particle on the real line is L2(R), the space 
of square-integrable functions on R in the Lebesgue measure. The operators corresponding to configuration and 
momentum variables act respectively as multiplication and differentiation operators. The Hilbert space is separable; 
an example of a countable basis are the harmonic oscillator eigenf unctions. 

To introduce polymer quantization on R, we start with the basis states 



The polymer Hilbert space is the Cauchy completion of the linear span of these basis states in the inner product 



where the quantity on the right-hand side is the Kronecker (rather than Dirac) delta. This space is clearly nonsepara- 
ble, and hence inequivalent to L2(R) [H!, [13, • Intuitively, building from the polymer basis states a single nonzero 
L2{M) state would require an uncountable superposition, and thus lead to an unnormalizable state in the polymer 
Hilbert space. Conversely, any state in the polymer Hilbert space has support on at most countably many points, 
and will thus represent the zero state in L2(M). 



II. POLYMER QUANTIZATION ON K 




(1) 



(2) 
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Next we define the action of the basic quantum operators. The position operator x acts by multiplication, 

{xi)){x) = XTpix), (3) 
and its domain contains the linear span of the basis states (JT]). The translation operators Ux, A G M, act by 

{Uxi^)(x) = i,{x + \), (4) 

and are clearly unitary. Formulas ([3]) and (jl]) are identical to those in L2{K)- In _L2(IR), the action of U\ is weakly 
continuous in A, and there exists a densely-defined self-adjoint momentum operator p such that p — ~'^\p\U\\y^_Q = 

—idx and l}\ = e^^^ . By contrast, in the polymer Hilbert space the action of Ux is not weakly continuous in A, and a 
basic momentum operator does not exist. 

The states in the polymer Hilbert space can be described as points in a certain compact space, the (Harald) Bohr 
compactification of the real line, and the operators introduced above can be described in terms of a representation 
of the Weyl algebra associated with the classical position and momentum variables [HI, [13, [i3| • There exists also a 
mirror-image quantization in which a a momentum operator and a family of translation operators in the momenta 
exist, but there is no basic position operator 14]. These mathematical structures will however not be used in the rest 
of the paper. 

As there is no basic momentum operator, any phase space function containing the classical momentum p, most 
importantly the Hamiltonian, has to be quantized in an indirect way. Following we fix a length scale A > and 
define 

P=^(^^A-?7r'): (5a) 



?^^{2-Ux-Ul). (5b) 
The Hamiltonian operator that corresponds to the classical Hamiltonian H — ^p^ + V{x) is then 



H^^i'i-Ux-Uh^V, (6) 

where V is assumed so regular that V can be defined by pointwise multiplication, {yip){x) = V{x)^{x). In ^2(1^), 
the A ^ limit in ([5]) would give the usual momentum and momentum-squared operators —idx and —9^, and the 
kinetic term in ([5]) would reduce to — ^9^. In the polymer Hilbert space the A ^ limit does not exist, and A is 
regarded as a fundamental length scale. For A ^ 1, one expects the polymer dynamics to be well approximated by 
the Schrodinger dynamics, and certain results to this effect are known |,[il,[li. 

Although the polymer Hilbert space is nonseparable, the dynamics generated hy H ^ breaks into separable 
superselection sectors. To discuss this, it is convenient to introduce a Dirac bra-ket notation in which the basis state 
-0^ © is denoted by |^). From ^ and dH) we then have 

(Ai|/i')-W, (7) 

*Im> = mIm), (8) 

;7a|m> - Im-A). (9) 

The action of H on gives a state with support at /i, /i — A and -t- A. The time evolution of thus has support 
only on the regular A-spaced lattice {/x-l-nA | n S Z}. This means that the time evolution breaks into superselection 
sectors, where each sector has support on a regular A-spaced lattice and is hence separable, and the sectors can be 
labelled by the lattice point /i € [0, A). The time evolution of any given initial state will consequently have support 
only on a countable union of A-spaced lattices. The upshot is that the time evolution of any separable subspace is 
restricted to a separable subspace. Thus, even though the polymer Hilbert space is nonseparable, the fundamental 
length scale A and the choice of an initial state or an initial separable subspace will result in quantum dynamics that 
takes place in a separable Hilbert space. 

We note that on a fixed A-spaced lattice, the Hamiltonian fG]) agrees with a conventional discretisation of 
Schrodinger's equation by the replacement 

^"{Xn) ^ {lpn+1 - 2VA> + i'n-l) ■ (10) 
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This suggests investigating versions of ([6]) in which the kinetic term is replaced by an operator that, in the finite 
difference approximation context, is higher-order accurate in A. The discussion of superselection sectors would extend 
to such versions with only minor changes. As the main interest in the present paper concerns singular potentials 
rather than higher order accurate discretizations of the second derivative, we shall work with 

In summary, the restriction of the polymer dynamics into any of its superselection sectors is mathematically equiva- 
lent to a conventional discrete approximation to the continuum Schrodinger equation on the corresponding equispaced 
lattice. The conceptual difference is, however, that in the polymer theory the lattice spacing is regarded as a new 
fundamental scale. 



III. THE RADIAL COULOMB PROBLEM ON A LATTICE 



1 



0) sector yields 



(11) 



Reduction of Schrodinger's equation with the Coulomb potential to the spherically symmetric {I 
the radial Hamiltonian operator 

H - 

acting in the Hilbert space of square integrable functions on (0, oo) in the measure dr [T7| . (For numerical convenience, 
the radial coordinate r has been chosen as twice the Rydberg radial coordinate.) J?rad has a one-parameter family 
of self-adjoint extensions, each characterized by a boundary condition at r — > [1^. The conventional choice for the 
self-adjoint extension is to assume the three-dimensional eigenfunctions in L2(lR'^,d^x) to be bounded at the origin: 
in terms of the rescaled wave functions on which -ffrad acts, this means that the wave function vanishes at r ^ p^ . 
The spectrum then consists of the positive continuum and the negative discrete eigenvalues 

en--4^, n=l,2,... . (12) 

We shall return to the other possible choices of the self-adjointness boundary condition in section FVl 

To introduce a polymer counterpart of -ffiad along the lines of ([5]), we need to address the positivity of r, the 
boundary condition at r = 0, and the singularity of the potential at r ^ 0. Suppressing the singularity issue for the 
moment, the first two issues can be solved by extending r to negative values: denoting this extended coordinate by 
x S K, and assuming that the discretized potential is symmetric under the reflection x ^ —x, we require the states 
to be antisymmetric under x i— *■ —x. In the two superselection sectors in which the lattice points are at respectively 
Xn — n\ and a;„ = (n + i)A, n £ Z, this antisymmetry condition just chooses the antisymmetric states. The 
remaining superselection sectors are pairwise coupled by the antisymmetry condition. However, since the potential 
is by assumption symmetric under x ^ —x, the energy eigenvalues can be found within each sector without using 
an antisymmetry condition, and the antisymmetric eigenstates are then obtained by just taking appropriate linear 
combinations. 

Let us return now to the singularity of the potential at a; — 0. Although this issue only arises in the single 
superselection sector that has a lattice point at a; = 0, we wish to give a prescription that handles all the superselection 
sectors in a unified manner. Let Xn denote the lattice points. Since 

sgn(a;) ^ ^ d{^/\x\) ^^^^ 
dx ' 



we can represent (sgn(a;)) /(yl^) by a finite-difference version of the derivative, 

Sgn(x„) 1 / r 1 r r\ 

r. r J [V\Xn+l \ - V Fn-l|J • 



(14) 



Taking the square, this leads to the lattice potential 

1 



I , si 

^ ( vK+il - VN'i-ilj 



(15) 



which is well defined even for a;„ = 0. The resulting Hamiltonian operator can be written in terms of the fundamental 
translation and multiplication operators as 

\_ /_ 1 / . „ , _ „ N 2 

A2 



1 



(16) 
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FIG. 1: The solid line is the Coulomb potential and the dashed line is the lattice regularized potential (|17p for A = 0.1 

and its action on the basis state is 

H\n) = ^ (2|m) - - A) - Im + A)) - 1 [V\F+M - ^\^^y l/i). (17) 

The potential term in the last form of (|16p could have been arrived at by considering a phase space version of the 
identity (jl3p . namely 

(18) 

where x and p are canonically conjugate variables. This would be the route that led to Thiemann's regularisation of 
inverse triad operators in LQG f5l|. In the present context, we have arrived at (|16p simply by representing a derivative 
by a finite difference on a lattice. A comparison of the Coulomb potential and its lattice regulated version is shown 
in Figure [T] for A = 0.1. There is a striking repulsive modification near a; = 0, a result which in its quantum gravity 
incarnation leads to a bounded curvature at the Planck scale [8| . 



IV. SPECTRUM 



Let us focus now on the superselection sector with the lattice points Xn = n\, A € Z. As this is the sector in which 
resolving the singularity at a; = is necessary, we expect it will provide the most interesting test of the resolution 
proposal pS]) . 

We look for energy eigenstates in the form '^^CnlnX), where the coefficients c„ are subject to the normalizability 

condition |cn|^ < oo and the antisymmetry condition c„ — —C-n- The energy eigenvalue equation with eigenvalue 
e reads 



HY,Cn\nX) =eY,Cn\nX). (19) 

n 11 
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It reduces to the recursion relation 



c„ (2 - A/„ - A^e) = Cn+i + c„_i, 



where 



Suppose from now on that e < 0. The asymptotic form of the recursion relation (|20p as n — > oo is 



c„(2 - A^e) = c„+i + Cu- 



lt has the linearly independent solutions 



(20) 



(21) 



(22) 



(23) 



It follows [19( that the exact recursion relation (|20|) has only one linearly independent solution that does not grow 
exponentially as n — > oo, and this solution has the asymptotic form (j23p with the lower sign, and is hence exponentially 
decreasing in n. We use this observation to set up a shooting method for a numerical computation of the eigenenergies. 
We start at an initial tiq 3> l/(— eA), in which regime the asymptotic recursion relation (j22p holds, compute Cuq-i 
using and then compute cq using the exact recurrence relation (|20p . Because of the antisymmetry condition 

c„ = —C-m the eigenenergies are those for which cq — 0. The accuracy of the method is monitored by increasing the 
value of TiQ until the results no longer change to the desired accuracy. 

When A — 0.1, we find that the lowest 14 eigenvalues are such that the quantity k — l/\/— 4e is a few per cent 
below the lowest 14 positive integers, being thus a good approximation to the continuum spectrum (|12p . For the 
higher eigenvalues the numerics becomes slow, and we do not have an estimate of when k starts to differ significantly 
from integers. 

For smaller A the numerics becomes slower but indicates convergence towards the continuum eigenvalues from below 
as A 0. Figures [2] and [3] show plots of co as a function of fc for A = 0.01: within the resolution of the plots, the three 
lowest roots are indistinguishable from k = 1, k = 2 and k = ?>. Figures H] and [5] show the lowest two eigenenergies as 
functions of A for 0.005 < A < 0.3. 



V. SYMMETRIC BOUNDARY CONDITION 



The results in section HVl show that the singularity resolution method (|16p gives good agreement with the continuum 
results for the antisymmetric boundary condition. While this is the boundary condition that arises from the conven- 
tional treatment of the three-dimensional Coulomb problem [l3| , we now wish to consider the resolution method for 
the — l/|x| potential on the full real line in its own right. We must then find also the energy eigenstates that are 
symmetric under x ^ —x. 

Staying in the superselection sector in which the lattice points are a;„ — nX, A G Z, the numerical algorithm of 
section HV] can be adapted to the symmetric boundary condition by starting again at uq 3> l/(— eA) but computing 
now ci and c_i. The eigenenergies are those for which ci = c_i. 

The numerics is now considerably slower than for the antisymmetric boundary condition. When A = 0.1, we find 
that there is a ground state at e « —4.94, well below the continuum hydrogen ground state, and the first 12 excited 
states are such that the quantity k = l/\/— 4e is approximately 0.2 above the lowest 12 positive integers. When A 
decreases to 10"*^, the lowest eigenvalue decreases and shows no evidence of converging to a limiting value as A ^ 0, 
whereas the higher eigenvalues appear slowly to converge from above to — l/(4n^), n = 1,2,.... Eigenenergies of the 
ground state and the first excited state for selected values of A are shown in Table HI 

The behaviour of the ground state in the limit A ^ has a counterpart in the continuum theory. In the continuum 
theory, the solutions to the eigenvalue differential equation -ffrad'0 = that do not vanish as r — > 0-|- have a 
logarithmic singularity there, and it is not possible to single out a boundary condition at r — > 0+ by a 'symmetric' 
extension to negative r. Instead, the self-adjointness boundary conditions at r — > 0+ can be parametrized by a length 
scale L e MU {oo}, such that the eigenenergies are the solutions to the transcendental equation 




7 




FIG. 2: The coefficient co as a function oi k = 1/V— 4e for 0.98 < k < 2.2, with A = 0.01. The zeroes are near k = 1 and k = 2. 




8 



-0.24 -1 
-0.25 - 


' ■ • . . 


-0.26 - 




-0.27 - 




-0.28 - 




-0.29 - 




-0.3- 




-0.31 - 




-0.32 - 


1 1 1 1 1 1 



0.0 0.05 0.1 0.15 0.2 0.25 0.3 

X 



FIG. 4: The lowest eigeneriergy as a function of A for 0.005 < A < 0.3. The vertical error bar of each point is 10 Convergence 
to € = —0.25 is apparent as A — > 0. 
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FIG. 5: The second-lowest eigeneriergy as a function of A for 0.005 < A < 0.3. The vertical error bar of each point is 5 x 10 
Convergence to e = —0.0625 is apparent as A — > 0. 
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TABLE I: The energy eigenvalue as a function of A for the ground state (eo) and the first excited state (ei) with the symmetric 
boundary condition. 



A -£0 -ei 

10^ IM 0153 

10"^ 14.8 0.181 

10"^ 32.3 0.196 

lO"* 58.5 0.207 

10"^ 93.9 0.214 



where 

G{z) -.^^{1 + z)~ln\z\-l/{2z) (25) 

and is the digamma function [l^. The usual boundary condition is obtained with L = 0. For aU the other values 
of L, the eigenenergies are shifted downwards: in the limit L ^ 0+, the lowest eigenenergy tends to — oo as — 
while the higher eigenenergies tend to the L = values from above, with corrections that are proportional to L [l^ . 
Our symmetric boundary condition on the lattice produces thus eigenenergies that appear at A to be in qualitative 
agreement with the continuum eigenenergies at i — *■ 0, although the rate of convergence on the lattice is significantly 
slower. 

The singularity resolution method appears therefore to be in qualitative agreement with the continuum theory also 
for the symmetric sector. The contrast between the symmetric sector and the antisymmetric sector shows, however, 
that boundary conditions on the quantum states at the singularity have a significant effect on the spectrum even after 
the singularity has been regularized. 



VI. DISCUSSION 



In this paper we have discussed polymer quantization in the — potential on the real line. We resolved the 
singularity of the potential at a; = by representing a derivative by its finitely-differenced lattice version, by a 
technique that mimicks the regularization of the inverse triad operator in LQG ,5] . Focusing on an equispaced lattice 
with a lattice point at x = 0, our numerical simulations indicated that the energy eigenvalues in the antisymmetric 
sector converge rapidly to those of the conventional continuum Coulomb problem as the lattice scale A goes to zero. 
This is not unexpected, since antisymmetry on the full real line corresponds to the conventional boundary condition at 
the origin in the spherically symmetric three-dimensional Coulomb problem. In contrast, for the symmetric sector we 
found that the ground state eigenvalue appears to decrease without bound as A approaches zero, while the eigenvalues 
of the excited states appear to approach the eigenvalues of the conventional Coulomb problem from above. The 
singularity resolution method in the symmetric sector thus yields dynamics that is qualitatively similar to that in 
the Coulomb problem with an M77.conventional choice of the self-adjointness boundary condition at the origin, with A 
corresponding to the length scale associated with this boundary condition. 

We view these results as evidence that the singularity resolution method yields physically reasonable results for the 
— l/|a;| potential on the real axis, whether one regards the finite difference equation simply as a discrete approximation 
to the Schrodinger dynamics in L2(K) or whether one regards the lattice scale A as a fundamental length within polymer 
quantization. We emphasize that the symmetric and antisymmetric sectors were found to have qualitatively different 
spectra, where the antisymmetric sector produces the conventional continuum eigenvalues in the limit of small A. This 
shows that even after the singularity in the potential has been regularized, boundary conditions that one may wish 
to impose at the locus of the singularity can have a significant effect on the spectrum. 

We note in passing that our singularity-resolution technique may also be of interest as a numerical technique in 
the context of pure Schrodinger quantisation, as an alternative to numerical techniques that invoke the asymptotic 
form of the continuum solution near the singularity [20j . To explore this suggestion, one would need to compare the 
convergence properties of our scheme, as evidenced in Figures |4] and [5l to the convergence properties of the matching 
scheme of 

Our results may be viewed as supporting Thiemann's regularization of the inverse triad operator in LQG [H]- 
Furthermore, they suggest that a boundary condition at the classical singularity may have a significant role also in 
the loop quantum gravity context, both when evolving through a spacelike singularity [2ll . [23 | and when setting 
boundary conditions at a timelike singularity 23] . 

There are however at least three subtleties in this respect, (i) We focused the numerical simulations on a regular 
lattice that has a lattice point at the origin. Will the situation remain similar also on irregular lattices, and is there a 



10 



systematic control on the singularity effects when the lattice is refined [HI, [T^, H^l ? In particular, would a significant 
symmetry-antisymmetry distinction emerge also on irregular lattices? (ii) We introduced the polymer quantization 
after first reducing the continuum Coulomb problem to the spherically symmetric sector. If one wanted to discuss 
polymer quantization corrections to the Coulomb energy levels from a phenomenological viewpoint, the polymer 
Hilbert space should presumably be introduced already at the level of three independent spatial dimensions, (iii) The 
polymer Hilbert space utilized was the genuine physical Hilbert space of the system, and the energy eigenstates were 
simply the normalizable solutions to the eigenvalue difference equation in the polymer inner product. In LQG, the 
polymer Hilbert space is only the kinematical Hilbert space, and further issues may emerge when the physical Hilbert 
space for solutions of the Hamiltonian constraint is introduced [l, [1^, . 
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